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a smooth transition we correct a result obtained by Grishchuk recently. We 
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the perturbed metric was not evolved correctly through the smooth reheating 
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(reheating) transition. 
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I. INTRODUCTION 



The discovery of anisotropies in the Cosmic Microwave Background Radiation (CMBR) 
IHII was a major advance for Cosmology. Presently, a series of experiments is being con- 
ducted to measure these anisotropies on smaller angular scales. It is likely that the CMBR 
anisotropies origin from primordial fluctuations of the metric and matter fields which have 
been amplified later on. Therefore, the outcomes of the ongoing and planned observations 
may have important consequences on the theoretical models explaining the origin and the 
evolution of cosmological perturbations, and more generally, on our ideas about the early 
Universe. 

Among these ideas is the theory of perturbations generated quantum- mechanically during 
an inflationary stage. In the present work our attention shall be restricted to inflationary 
models with one scalar field. In the subsequent evolution of the Universe matter shall 
be described by a perfect fluid. With the help of various models definite predictions for 
the amplitudes and spectra of the different types of perturbations (density perturbations, 
rotational perturbations and gravitational waves) can be made. In particular, allowing 
perfect fluids only, one can show that rotational perturbations have to decay. Moreover, 
single scalar field inflation cannot seed rotational perturbations. 

The power spectrum of density perturbations in an inflationary Universe was first calcu- 
lated by Mukhanov and Chibisov ||. This first computation was confirmed in Refs. || 0. 
Essentially, the 'standard result' lies in the following: the closer the spectrum is to the flat 
(Harrisson-Zeldovich) spectrum, the larger the amplitude of density perturbations is in com- 
parison with the amplitude of gravitational waves |7]J| . On the other hand, if the spectrum 
is tilted away from scale-invariance, the contribution due to gravitational waves can become 
important. (For some specific examples of this behavior see, e.g. 0.) 

This result was recently challenged by Grishchuk in Ref. (K| who found a similar am- 
plitude for density perturbations and gravitational waves. However, his calculations were 



criticized by Deruelle and Mukhanov ||TT[. They argued that Grishchuk did not properly 
consider the joining conditions at the transitions 'inflation-radiation' and 'radiation-matter'. 
Then, Grishchuk published a comment |]T2[ in which he stated that the equation expressing 
the 'standard result' is mathematically wrong. Short after, this claim was contested by 
Caldwell [|I3| who expressed his agreement with Deruelle and Mukhanov. Finally, in the 
appendix of Ref. |14j], Grishchuk re-stated his criticisms of the 'standard result' in more 
details. 

The aim of this paper is to clear up this controversy and to study the influence of cosmo- 
logical transitions undergone by the Universe during its history on superhorizon cosmological 
perturbations. 

This article is organized as follows: the second section is devoted to density perturbations 
and gravitational waves in the gauge-invariant formalism. This section can be skipped by 
specialists. The third section deals with the synchronous-gauge formalism. It was claimed 
by Grishchuk that the synchronous gauge and gauge-invariant results differ by an arbitrary 
constant of integration. Thus, a systematic comparison of the two formalisms is made and 
we show that they are equivalent. In the fourth section, we study the so-called 'conservation 
law' for superhorizon modes. We show that the 'conserved quantity' ( often used in the 



literature is actually not empty as it was claimed in Refs. [12]. In the fifth section, we 
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turn to the question of the matching conditions. In particular, we re-derive the Deruelle- 
Mukhanov junction conditions with a different method. The sixth section is devoted to the 



study of the smooth transition of Ref. | ID| . We will argue that the joining conditions are not 
the essential point, but that the evolution of the metric perturbations through the smooth 
reheating transition was not done correctly. We correct the result, which is now the same as 
for a sharp transition. In the last section, we briefly present our conclusions, which confirm 
the 'standard result'. Finally, an appendix reviews how the initial conditions are fixed 
when the perturbations are quantum-mechanically generated. For the readers convenience 
we summarize the notation of the gauge-invariant formulation and the synchronous-gauge 
formulation in two tables at the end of the paper. 



II. GAUGE-INVARIANT SCALAR AND TENSOR PERTURBATIONS 



We will use the notation of [IS] except the signature of the metric and the harmonic de- 
composition of the gauge-invariant potentials. The line element for the Friedmann-Lemaitre- 
Robert son- Walker (FLRW) background plus scalar perturbations reads (c = 1) 

ds 2 = a(r]) 2 {-(l + 20)dr? 2 + 2B\ j &x i dr) + [(1 - 2^)7y + 2E\ Aj ]dx i dx j } . (2.1) 

The conformal time rj is related to the cosmic time t by dt = a(r])dr]. A dot denotes a 
derivative with respect to t, whereas a prime stands for a derivative with respect to 77. The 
tensor 7^ is the metric of the three-dimensional space-like hypersurfaces and the derivative 
I in Q2.1| ) is covariant with respect to 7^. The curvature of these sections is given by K and 
takes the values 0, ±1. 

One can generate fictitious perturbations by performing an infinitesimal change of co- 
ordinates which preserves the scalar form of equation ( |2.1| ) (let us note that this is not, as 
written in [I3~| on page 2439, a "general" infinitesimal transformation of coordinates): 



r ] = r ] + e(v,x k ), x i = ^ + 7%(7 7) x fc ) . (2.2) 

Two independent gauge-invariant scalar metric potentials may be constructed from the 
metric. Following Ref. |TSJ we take them to be: 

$Q = </> + ~[{B-E')a]', #Q = V- -(£-#')■ ( 2 - 3 ) 
a a 

In these expressions, we have chosen to extract the scalar harmonic Q(x l ) of the gauge- 
invariant variables from the very beginning. The function Q{x % ) satisfies the Helmholtz 
equation 

AQ = -k 2 Q, (2.4) 

where AQ = "f^Quti an d k is the comoving wave number. For convenience the quantity 
7i = a' /a is defined, which is related to the Hubble parameter H = a/ a = Jija. The 
perturbed Einstein equations can be expressed in terms of $ and alone (see below). 



Working in the gauge-invariant formulation of Bardeen [16] is equivalent to the longitudinal 
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gauge (B = E = 0), which fixes the constant time hypersurfaces to be the hypersurfaces 
with vanishing shear. 

In most inflationary models matter is described by a scalar field ip = foiv) + 8ip(r],x l ). 
The background energy density and pressure are 

po = ^ + v%o), P° = l|r-^ o) ' (2 ' 5) 

where V((p) is the potential of the scalar field. If (p' 7^ 0, the covariant conservation of the 
energy-momentum tensor provides the Klein-Gordon equation 

<p% + 2H<p' + a 2 V tip = . (2.6) 

The gauge-invariant scalar field perturbation is Sip^Q = 5ip + ^p' (B — E'). The linearized 
Einstein equations for a fixed mode k may be written in terms of gauge-invariant quantities 
only: 

-3H{H<S> + - k 2 ^ + = |[-(v^o) 2 $ + ^o(^ (gi) )' + aV»>^ (gi) ] , (2.7) 
H<5> + = -<^o^ (gi) , (2.8) 
$ - * = , (2.9) 

(2V! + H 2 )§ + H& + *" + 2W - - ifc 2 ($ - *) = 

3 

|[-(^) 2 $ + ^(^Wy - a 2 VJ^\ , (2.10) 

where k = 87rG. Everything can be expressed in terms of a single gauge-invariant quantity 
since equation ( |2.9| ) tells us that $ = \I>. 

If <p' = the background solution is the de Sitter spacetime (p = —Po)- In that case, 



the solution of the system fl2.7f) - ( p . 1 Q|) is $ = 0: there are no scalar metric perturbations. 
This does not mean that the scalar field cannot fluctuate but that these fluctuations do not 
couple to gravity. 

For (p' 7^ the linearized Einstein equations (|2.7| ) - Q2.10| ) reduce with help of the 
Friedmann equation and the Klein-Gordon equation Q2.6Q to the equation |l5j 



$" + 2{H - ^f)$' + [k 2 + 2(H' - - 2iT)]$ = . (2.11) 
For K = 0, the introduction of the new variables[] 



allows us to express equation (|2.11|) in the form: 



1 The perturbation a is related to u of Ref. [15] by a = (/•t/3) 1 / 2 u. 
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a" + (k 2 - ^-)o = . (2.13) 
9 

For modes k 2 <C 9" j 9 the solution of Eq. (|2.13 ) may be expanded in powers of k 2 . At 
the leading orders we obtain 



a = A 1 (k)9 J-(l-k 2 f 9 2 f -dfjdf) + 0(k A )^ d V + (2.14) 
A 2 (k)9 (l-k 2 J^J V 9 2 df}drj + 0(k 4 



Since 9 — > oo for a — > in general, is the arbitrary constant in front of the regular 
(growing) mode and A 2 a constant associated with the singular (decaying) mode?]. 
For gravitational waves the line element reads: 

ds 2 = a(ri) 2 {-dr] 2 + [ 7ij + h^dx^} , (2.15) 

the tensor hij being symmetric, traceless and transverse. The tensor sector is gauge- 
invariant. We write = h^Qij, where Qij is a symmetric, transverse, and traceless spher- 
ical harmonic, and h gw is the amplitude of a gravitational wave. The decay of the amplitude 
due to the expansion of the Universe is taken into account by defining /i gw = 0(77)^^(77). 
The equation of motion for /i^ reads: 

AC + (k 2 - ^>gw = . (2.16) 

A physical interpretation of Eqs. ( |2.13j) and (|2.16j ) is parametric amplification of the per- 
turbations while the Universe is expanding JT7J. The scale factor plays the role of a 'pump 
field' and the 'interaction' between the background and the perturbations is described by 
the 'potentials' 9"/9 and a" /a. The potential of the scalar perturbations involves not only 
the scale factor but also the derivatives of a(rf) up to 0^(77). 

For simple models where the scale factor is given by 0(77) = lo\i]\ 1+f3 , the exact solution 
to Eq. ( |2.16| ) can be found. It reads: 

/V = (^)V^r(A0/ W2 (&77) + A*T(k) J^ +im (kr,)}, (2.17) 

where J±(/3+i/2) are Bessel functions. 

For an arbitrary scale factor, as for the scalar metric perturbations, the solution for 
modes k 2 a" /a is given for gravitational waves by replacing 9 with a in Eq. ( |2.14[ ): 

/i gw = Af{k)a + Af(k)a + 0(k 2 ). (2.18) 

Now Af w corresponds to the regular and v4f w to the singular mode. It should be noticed 
that the v4 gw 's differ from the A gw 's, although they are connected. At superhorizon scales 
the dominant mode is constant in time, independent of the matter content of the Universe. 



2 An example where both modes are singular as a — > is provided by the model with scale factor 
behavior a oc for —2 < j3 < —1 and 77 < 0. However, this inflationary model violates the 

weak energy condition pq + po > and cannot be realized with a single real scalar field. 
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III. SCALAR PERTURBATIONS IN SYNCHRONOUS GAUGE 



Let us describe the scalar perturbations using the class of the synchronous gauges. In 
order to make contact with previous works, we will use the notation of Ref. |I0| . There the 
quantity H is denoted by a, a = TL = a' /a and the function 7(77) is defined by: 

7 = 1--^- 3.1 

This function reduces to a constant for scale factors which are proportional to a power of the 
conformal time. The function 7 is zero if (p' = (de Sitter). In Ref. [|Kj, the line element 
is written as: 

ds 2 = a 2 (r]){-dr] 2 + [(1 + fcQ) 7y + ^-L^Q^.JdjrW}. (3.2) 

Choosing the synchronous gauge (SG) means setting the perturbed lapse and shift functions 
to zero. The notation of |ID| is related to the notation of by 

= 0, S = 0, V = ~0, E= 2(k* h -K) Q - (3 ' 3) 
Thus, the gauge-invariant variables expressed in the synchronous gauge are given by: 

¥SG) = -V h -W^K h '^ (3 ' 4) 

^ SG) = ~2W^K) {h ' l ' + ah[h (3 ' 5) 
StpWW = fa - 2(jfe^)«) , (3-6) 

where the scalar field is written as: tp = ^o(^) + fi(v)Q- Inserting these formulas into the 
system ( |2.7| ) - ( |2.10| ) provides the correct perturbed Einstein equations in the synchronous 
gauge (see Ref. flQfl), namely: 

h" + 2ah[ - (k 2 - K)h = , (3.7) 

~ h ' + k 2 - K kl = K(f ' oifl ' 
3ati - ah\ + (k 2 - AK)h = re(^i + c?ipiV )V ) , (3.9) 
- h" - 2ati + Kh = K(^oVi - a 2 fiV v ) ■ (3.10) 

This shows that, as expected, the two formalisms are completely equivalent. It is worth 
noticing that equation ( |3.7| ) expresses the fact that \I/ = $. This is due to the vanishing 



longitudinal (anisotropic) pressure (see Ref. |L0 



Let us now consider the question of the residual gauge. It has been known for a long time 
that the condition h 0fl = does not fix the gauge completely. This condition is preserved 
under the change of coordinates 

' = §«"/ T-f« K - (3 - n) 
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where C and D are arbitrary constants for a fixed mode k. The corresponding changes for 
h and hi are given by: 



h = h + C-, hi = h + k 2 C f-^- + k 2 D . 
a J a 



(3.12) 



According to Ref. flOA , one can thus introduce two 'residual-gauge- invariant' quantities u 
and v defined by: 



u = h' + a^h, 



k 2 

v = h', h . 

a 



(3.13) 



The relations between the gauge-invariant quantities and the residual-gauge-invariant quan- 
tities can be expressed as: 



(3.14) 
(3.15) 



In Ref. [10] it has been claimed that u and v are genuine gauge-invariant quantities [i.e., 



invariant under the transformation (|2.2f) l. This claim is not correct. Indeed, a direct check 
shows that: 



u = u + 2a f + 2a 2 £° 

v=v+2k 2 (e-e) . 



(3.16) 
(3.17) 



The equations (|3.14j) - ( |3.15[) just give the value \I/ and <3> calculated in the synchronous 
gauge, ty( SG \ $( S ' G ). It does not come as a surprise that what remains in a fixed gauge from 
the gauge invariance is simply the residual gauge invariance. 

In Refs. [0] and 0], it has been argued that equation ( [2.1 1| ) is incorrect since it could 
be expressed as a combination of the derivatives of the correct equations and hence would 
contain a non-physical constant. This claim is incorrect as well. Let us demonstrate why. 
If we insert the expression of $ in terms of u and v [equation ( |3.15| )1 in equation ( 2.11| ), we 
find that this one transforms identically into: 



2k 2 a 



[k 2 u + a(v' + 2av)]" + ^(2^-+ ' 

k ™ A 



a' 
lx 2 



- -)[k 2 u + a(v' + 2av)}' 
2ct7 2 



+ 



a 



2a 



a'i 



+ 



a 



2 L 



a. 



av 



7 



(- + -) 



u 







a 
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(3.18) 



Of course we must also take into account the equation \l/ = $ which is, in terms of u and v, 
equal to: — k 2 u = a(v' + 2av). Therefore, we see that from the gauge-invariant formalism, 
we can reduce the whole problem to a set of two coupled first order differential equations 
(of course, this could have been done directly in the synchronous gauge) for the variables u 
and v : 
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— k 2 u = a(v' + 2av) 



u 



v = i(—y. 



(3.19) 
(3.20) 



From this system, we can generate two decoupled second order differential equations: 

, (3.21) 
. (3.22) 



u" + u'(2a^ - -) + u[k 2 - 2a' -a 1 -- /7 ' V1 



7' 



7 



7 



v" + u'(2a - — ) + k 2 v + v(2a' - 2a- 7 



7 



7 



The last equation could have been guessed from the very beginning by inserting the ex- 
pression & SG *> = \&( SG ) = (cry) / '(2k 2 ) in equation ( |2.11| ). The variable change (see Ref. 
ITT 



«\/7 

U = jj, 



(3.23) 



in formula (|3.21|) allows us to obtain the correct equation for //: 



awr/ 



(3.24) 



Therefore, we have proved that the gauge-invariant framework leads to the same (correct) 
equation of motion for the variable \x as the calculation in synchronous gauge. Let us note 
that the residual-gauge-invariant variable /i is nothing but the value of the gauge-invariant 
variable 



v M = a(6<pW + ^$ ) 
ri 



(3.25) 



which was defined by Mukhanov in Refs. [18,15], expressed in the synchronous gauge. In- 
serting the synchronous gauge values of Scp^ and $ into (|3 .25 ) and using ( |3.19| ) yields: 



(SG) 



'2k 



(3.26) 



In order to be as complete as possible we examine what was the mistake of Refs. [|T2| , [14 
From the equation ( |3.14j ) and the definition of /i, we can re-write equation ( |3.20| ) as: 



$ (SG) = 



cry fi 



2k 2 a^' 



(3.27) 



If one inserts this expression in formula ( 2.11]) , we obtain a third order differential equation 
for /x: 



2 ( a V^)"i v ( a V / 7)'r..// , r/„2 ,( ' 



-w 



a^7 a^7 
which can be integrated and leads to: 



a V7 



(3.28) 
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H" + [k 2 - 



(3.29) 



where X is a constant of integration. Comparing equations Q3.24| ) and ( |3.29| ) shows that 
X = 0. 

Let us demonstrate why the introduction of \i via fl3.27| ) and equation ( J2.11[ ) seems to 
give X ^ at first sight. Consider a set of two first order differential equations [they play 
the role of equations (|3.19|) and (|3.20|) , the latter is equivalent to (|3.27|) 1: x' = —y, y' = x. 
From them we can generate two second order differential equations [they play the role of 
equations (|3.21|) and (|3.22|) , which is nothing but equation ( j2.11| )1: y" + y = and x" + x = 0. 
But we can also insert y' = x into the second of the two last equations [as we inserted ( |3.27| ) 
into ( p.ll| )1. We will obtain a third order differential equation, y'" + y' — 0, which may be 
integrated to yield y" + y = X . There is no harm to do that as long as we do not forget to 
use equation y" + y = and therefore X must be equal to zero. 

However, in Refs. fl2,14] the derivation of equation ( 3.27 ) was not given, instead it was 
assumed that ( |3.27|) may be considered to be the definition of fi. In that case, this 'new' 
\i has nothing to do with the [i defined previously and has not to satisfy Eq. (|3.24j ). The 
equation of motion for this 'new' variable /i is Eq. ( |3.29| ). This definition of fi is unique 
up to a shift fi — > // + Ya^py only, where Y is an arbitrary constant for a fixed mode k. 
Demanding that \i should fulfill ( ggg ) fixes Y = X/k 2 . The choice Y = 0,X ^ that 
implicitly was made in |T2| , ^4| is inconsistent. 

In conclusion, let us emphasize again the main result of this section. The gauge-invariant 
formalism and the synchronous gauge formalism are completely equivalent for all values of 
k including the zero-mode. Using one or the other (or any further gauge) is only a question 
of taste or of prejudices. 



IV. THE CONSTANCY OF C FOR SUPERHORIZON MODES 

In this section we turn to the study of the so-called 'conservation law'. Let us start by 
introducing matter that can be described by a hydrodynamical equation of state, in order 
to follow the evolution of the perturbations from reheating till today. 



A. Perfect fluids and the 'standard result' 



Assume that anisotropic stresses can be neglected, thus matter is described by a perfect 
fluid. The pressure of the perfect fluid is related to its energy density by the equation of 
state p = p(p, S). S is the entropy per baryonQ, n is the density of baryons. For a reversible 
expansion of the background (there are no unbalanced creation/annihilation processes) the 
entropy per baryon is constant in time. Due to the isotropy of the background, the first law of 
thermodynamics reads d(p/n) = —pd(l/n), thus TdS = from the second law for reversible 



Of course, this makes sense when baryon number is conserved only. 
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processes. Therefore, the background equation of state has to be isentropic (VjS 1 = 0), thus 
we may write p = p(p) = w(p)p for an adiabatically expanding background. 

The pressure perturbation reads 5p = c 2 5p + t5S, where c 2 = (dp/dp) s is the isentropic 
sound speed and r = (dp/dS) p . For an adiabatically expanding background p' = c 2 s p' . 
Below we use the relation c 2 = p'/p', thus the following results hold true for negligible 
entropy production only. 

The equations of motion for the gauge-invariant metric potentials for a perfect fluid are 
(see Ref. |T5fl) $ = * and 



$" + 3(1 + c 2 s )H& + [2H' + (1 + 3c 2 s )(H 2 - K)]$ + c 2 s k 2 <& = ^a 2 r5S , (4.1) 

except for the de Sitter universe, where $ = 0. Written in terms of the variable a, which is 
defined by (|2.12|) , the latter equation readsQ: 



o" + [c 2 k 2 - 6 -pj a = —a^rSS . (4.2) 

For isentropic perturbations the leading order solution is easily obtained to be a = 
C\9 J drj/8 2 + C2O, whereas the next to leading terms differ from the solution ( |2.14j ), be- 
cause c 2 may be time dependent. [There should be a factor c 2 in the fj integration in (|2.14j ).l 
If w and the sound speed are constant an approximate solution of a is not of much use, 
because then the exact solution to all orders in k can be given in terms of Bessel functions. 
Two examples of a perfect fluid are the radiation fluid (w = c 2 = 1/3) and dust (w = c 2 = 0). 
These examples have vanishing entropy perturbations (5S = 0). In general, 5S does not 
vanish for more than one fluid. The scalar field if is another form of matter that can be 
described by a perfect fluid, see (|2.5|) . Formally, the sound speed is defined as above. With 
help of the Klein-Gordon equation (|2.6p it reads 



= -K 1 + H) ■ (4 - 3) 

Now we define the 'entropy perturbation' through t5S = 5p — c 2 5p. With the expressions 
for dp and Sp and using Eqs. ( p.7|) and (|2.8| ) we arrive at 

^a 2 r6S = (l - cjfa,)) (3JT - A; 2 ) $ . (4.4) 

We obtain the scalar-field equation of motion for the metric potential, Eq. ( |2.11| ), by inserting 
( |4.4|) into (|4.1|) and replacing c 2 by (|4.3|) . Thus, we may study scalar metric perturbations 
by means of Eq. ( |4.1| ) from their generation during the inflation epoch to its observation in 
the CMBR today. However, during reheating [plj the scalar field may oscillate, thus ip' has 
zeros and Eq. ( |2.11| ) is singular at these points. This situation has been discussed in 
Let us assume in this work that this does not happen. 



4 With the definition 6 = l/a[p /(p +p )] 1/2 (l - SK/npoa 2 ) 1 / 2 Eq. (gj) holds true for all values 
of K. 
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We are now able to derive the 'standard result', i.e. the amplification of $ during the 
reheating transition. From comparison of Q2.14 ) with the solution of (4.2) it is clear that 



the leading superhorizon term of the solution does not depend on the sound speed c s . After 
some time the decaying mode is unimportant and the leading superhorizon term, using 
( P13 ), reads: 



$ - 2 lk) ^ J a (1 + w)dr? = Al(A;) ^ y a y-~ w) dn • (4 - 5) 

For a power law behavior of the scale factor, i.e. a oc |^| 1+/3 , the equation of state is given 
by w = (1 — /3)/[3(l + /?)]. Then, the 'growing' mode is constant in time: 

®^l Al{k) WV3 {1 + wh (4 ' 6) 

Let us assume that the evolution of the scale factor may be described by such a power law far 
away from the transitions inflation-radiation and radiation-matter. In between it changes 
smoothly. The ratio of the values of $ during inflation and matter domination (again far 
away from the transitions) is then given by 

$™ _ 2 2A + 3 1 ^ 2 1 

$i ~ 5 A + 1 1 + Wi ~ 5 1 + ' 1 ' j 

where /3 m = 1, ~ —2, and u>j ~ —1. Therefore, scalar perturbations are magnified by a big 
factor during the reheating transition. For de Sitter spacetime the amplification coefficient 
goes to infinity. This simply expresses the trivial fact that $ goes from zero to a constant 
without prejudice to the numerical value of this constant. Note that we can not conclude 
from this argument that $ m is large because $j — > as Wi — > — 1! The real (absolute) value 
of $ after the transition can be only known after having determined the initial conditions 
from the quantization of density fluctuations, that is to say after having fixed A\. 



B. Definitions and use of ( 

As was first recognized by Bardeen, Steinhardt and Turner and further elaborated in 
the equation of motion ( (4.1| ) has a first integral for isentropic modes (5S = 0) that 
are much larger than the Hubble scale, i.e. fc phys = k/a <^ H. Following Ref. JL5] we define 



2?T 1 $ / 
CS 3 



1 



w 



(4.8) 



which was introduced by Lyth [^] originally (a quantity differing by terms Oik 2 jH?) only 
was used by Brandenberger and Kahn |22|). £ essentially is the perturbation of the intrinsic 
curvature in the comoving gauge [2~3l|. Its first derivative reads 



Lc = 2 1 

H 3(l + w) 



K /l + 3w 



n 2 



2H 



$' + 3(c 2 



, . k a 2 r5S 



(4.9) 



where we used the equation of motion (|4.1| ) and the equations 
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w 



3W(1 + w)(w - c 2 s ) 



Ti' 



1 + 3w 



{Ti 2 + K) 



(4.10) 



Let us note that (4.9) is an equation for a fixed mode k, i.e. the large scale limit k/Ti — > 
in this equations means to make 7^ large. Thus, C is constant in time for superhorizon 
modes k/Ti, <C 1 if and only (i) if = 0, (ii) entropy perturbations are negligible, and (iii) the 
perturbation is given by the regular mode only. The conditions (i) and (ii) are obvious from 
( fl.9| ). The last condition means that the last term on the r.h.s of ( |4.9| ) vanishes in the limit 
k/Ti — > 0. The decaying mode is $ oc Tia/(a 2 9) oc Ti/a 2 at leading order. Thus, the r.h.s 
is proportional to (k/Ti) 2 Ti/a 2 , which blows up in the limit a — > 0, except for very special 
casesf]. Therefore, in general, one has to exclude the decaying mode in order to make use of 
the 'conservation law' C' = 0. In fact C is nothing more than a first integral of the equation 
of motion fl4.1|) when the conditions (i) - (iii) are fulfilled. 

The above definition of £ differs from the original definition |]6|j24 







BST 



Sp 



3 Po + Po 



(4.11) 



Cbst is a hypersurface-independent quantity |24 
metric potential it reads: 



Written in terms of the gauge-invariant 



Ti 2 



,BST 



3 (1 + w)(Ti 2 + K) 



! . , K 1 f k' 



- $ 



(4.12) 



From its time derivative 



1 



o ; Cbst 



7i 2 



3(l + w)(H 2 + if) 



/-I 



^ . _ Ka 2 r5S 

$ + $) H 

J 2 Ti 2 



(4.13) 



the constancy of Cbst follows if and only (i) there are no entropy perturbations, and (ii) 
there is the regular mode only. 

Besides the advantage of Cbst over ( to be conserved even if K ^ 0, Cbst is a hypersurface- 
independent measure of the metric perturbations including regular and singular modes (sin- 
gular in the limit k/Ti — ► 0), whereas C does not measure singular modes, if one takes the 
leading order contribution into account only. This can be easily seen by rewriting C in terms 
of o: Whatever the coefficient C 2 in front of the decaying mode is, it does not enter into C 
in the leading order in k/Ti, because 



C = 9 2 



(4.14) 



for K = 0. On the other hand 



5 In the model given by the scale factor a oc r) 1+ @ with —1 < (3 < —1/2 and r\ > the decaying 
mode does no harm, because (k/H) 2 Ti/a 2 oc r]~ 2 ^~ x . However, these models lead to equations of 
state where 3j>o > Po, which does not give rise to inflationary expansion. 
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C: 



BST 



a — 



e 



3H 



(J 



(4.15) 



(for any value of K) does depend on both C\ and Ci- Let us in the following discuss the 
properties of C(&)- 

In Ref. |12[] it has been claimed that "the conservation law ((ti) = C{tf) degenerates to 
an empty statement = 0". In order to explain the line of reasoning of |12| and to show 
where the argument fails, let us restrict the discussion to matter in form of a scalar field. 
For an isentropic perfect fluid the line of reasoning is analogous. Since ( is a gauge-invariant 
variable its value is the same in all gauges. Let us compute it in the synchronous gauge. 
Inserting Eq. (jjEZTp in Eq. (|4|) we find that: 



C 



a 2 7(- 



2fc 2 a 2 7 ay/ 7 ) 



(4.16) 



In Ref. Eq. (ggg ), re- written as 



a 2 7 



[a 2 7 (. 



// 



<V7 



ttyff 



X, 



(4.17) 



was inserted into ( 4.16 ) to obtain ( ~ X/(2k 2 ) in the limit k — > 0. According to Ref. ||12|| . 
there is no mean to know that X = in the gauge-invariant formalism. This would explain 
( = const. / in the considered limit. On the other hand, always according to Ref. |12j| , the 
synchronous gauge formalism could tell us that X = 0. Thus, we would be 'betrayed' by the 
gauge-invariant formalism in which X would appear. Only computations in the synchronous 
gauge formalism could reveal that £ = X/(2k 2 ) = 0. In Refs. [|T^,0], the fact that ( = at 
the 'leading order' k~ 2 was misused as a proof that it is not possible to use the quantity Q 
to learn something about the behavior of density perturbations. 

We have shown in Sec. Ill that the synchronous gauge formalism and the gauge-invariant 
formalism are completely equivalent and that X = in both approaches. Thus, there is no 
risk of confusion at all computing ( in the gauge-invariant formalism. 

According to Ref. |13j (p. 6), the conclusion of Ref. jrj] occurs because "the k — > limit 
has not been taken consistently" . Let us show that this is not the case. Following Ref. [12 



we assume X / because we now define \x by Eq. ( |3.27|) , but remember that this definition 
of /x is not unique. In that case the equation of motion is given by Eq. (|3.29|) . Let us expand 
the solution \i{r\) in powers of k 2 : 



ttyf) 



A^k) [l-k 



+A 2 {k) 



+X 




1 - k 



a 2 -/df]dr] + 0{k 4 ) 
rv 1 



f ^ f + dr] + ( 418 ) 

l-k 2 f±- J V a^df/dfj + 0(k A )\ dfjdr] , 



where Ai(k), ^(/c) are arbitrary integration constants fixed by the initial conditions. If we 
compare Eq. ([4.16|) with Eq. ([4.17|) we deduce that: 
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2k 2 2a^y 
and we find the beginning of the series giving (: 

X lf 1 - r drj r 1 rv 



C = 4"^. (4-19) 



C =2F-2^« + * ( ' ;) /^ + A 7^/"^ d '' d " + 0(P) ) ' < 42 °' 



This confirms that the limit had been taken properly in Ref. [12] and that the first term in 
the series contains X/k 2 . This result does not imply ( — > oo as k — > 0, because X = X(k) 
and we know nothing about the k dependence of X. However, there is no argument telling 
us that Ai(k) or A 2 {k) are subleading compared to X{k)/k 2 . 

Taking into account the correct value X = 0, Eq. ( |4.20| ) reads: 

c = 4-T/r^ + °<' ;2 >- ("I) 

2 2 J {dy/j) 2 

The same expression could have been obtained directly within the gauge-invariant formalism. 
It is necessary to push the expansion in k 2 one step further because the leading order term 
of the decaying mode does not contribute in ( [4.14j ). Inserting Eq. fl2.14|) in Eq. Q4.14|) , we 
obtain: 

( = A 1 - A 2 k 2 f dr]6 2 + 0{k 2 ) . (4.22) 



Comparison with Q4.21| ) yields A^k) = —2A\{k) and A 2 (k) = 2k 2 A 2 {k). We would have 



obtained the same relations between the initial conditions by inserting ( |4.18|) into ( 3.27) and 



comparing the result with ( |2.14 ). From ( |4.22| ) it is seen that ( does not vanish and is not 



constant in general. The second term of ( |4. 22| ) is the decaying mode. We emphasize that it 
is crucial to neglect the decaying mode if one makes use of 

C ~ At . (4.23) 

Let us show how to make use of the constancy of £. Assume conditions (i) - (iii) are 
fulfilled. An example is the 'standard' scenario: The initial perturbations are provided by 
quantum-fluctuations during inflation. At the first horizon crossing we denote them ^(tj). 
Let k be a mode which re-enters the Hubble horizon after equality between matter and radi- 
ation. The theory of quantum fluctuations shows that <&k{ti) is due to isentropic (adiabatic) 
fluctuations. The decaying mode is negligible short after the first horizon crossing. Entropy 
perturbations due to the reheating transition and due to the transition from radiation to 
matter affect much smaller scales. We know that the growing mode in $ is a constant on 
superhorizon scales if the scale factor obeys a simple power law behavior. In the inflationary 
stage a (xt p where p 3> 1, in the matter stage a oc t 3 ^ 2 . Thus, with w(t m ) m 

^ f^^Mfi) - l^itm) • (4.24) 
1 + w[ti) 6 

During inflation w ~ — 1 and therefore the large amplification 
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$ fc (0 ~ |c * I n ■ Sn ^ftO ( 4 - 25 ) 
5 5 1 + iy(tjj 



follows, which is the same as 

To finish this section let us consider the concrete model studied in Ref. WM (see the 
appendix of that paper). It consists in a transition from one power law scale factor, a(t) = 
ait Pl for t < ti, to another power law scale factor, a = a 2 (t — t*) P2 for t > t\. At the 
transition a and H are continuous, whereas H jumps. Assuming the constancy of C = Co we 
may integrate the definition of C fl4.8|) in terms of cosmic time to obtain the evolution of the 
Bardeen potential 



For t < ti we find: 



S(t) = -Co- / ^ + -C (4.26) 
a Ju ti a 



$(t < ti) = -^[1 - (-) 1+Pl ] + — ^— C. (4.27) 



Consistent use of the constancy of £ = Co requires the vanishing of the decaying mode in 
Eq. (]427D and therefore C = Coait- +pl /(pi + pf). We obtain: 

$(t < tl ) = T ^-. (4.28) 

For the de Sitter case ( l p\ = oo') we recover that $ = 0. Let us compare this result 
with the one of Ref. ||14|| . We do not agree that "The initial value of the potential is 
®(ti) — Co + H(ti)C/a(ti), ..." and that "The constant C could be set to zero from the 
very beginning." as it is stated on page 31. This would lead to $(tj) = Co an d, would imply 
a non-vanishing Bardeen's potential for the de Sitter spacetime. Instead the initial value is 
®(U) = Co/(l +Pi) which is clear from Eq. ( |4.28| ). 
For t > ti, we obtain the solution: 

$(t > h) = _^_fPly+P*(Jl_)i+Ps + _ jLji+^i (4 29 ) 

v ' 1+piPi y t-tJ l+p 2 [ K t-tJ 1 v ; 

The Bardeen variable is continuous at t — t\ (since we have integrated a Heaviside function). 
Long after the transition, the previous relation reduces to: 

$(t>f 1 )~_^. (4.30) 

1+P2 



Therefore, we reach the conclusion that the amplification coefficient is given by: 

${t > ti) _ 1 +pi 
$(t < ti) ~ l+p 2 



(4.31) 



Formula ( [4.31D should be compared with the third equation after Eq. (87) of Ref. [I3|. In 



this expression the missing term 1 + pi in the numerator is due to the incorrect assumption 

m) = Co- 



lt 



V. A SHARP TRANSITION JOINING CONDITIONS 



There are two physical situations where a sharp transition in the equation of state is 
a good approximation. The first one is the study of the behavior of superhorizon modes 
(fc phys <C H). The duration T of processes like the change from the radiation dominated to 
the matter dominated universe, reheating at the end of inflation, or recombination typically 
take several expansion times, i.e. T ~ H^ 1 . Superhorizon modes change on much larger 
time scales, i.e. l/k phys , thus they see a sharp transition. These sharp transitions may 
violate the second law of thermodynamics, e.g. at equality radiation entropy is destroyed 
instantaneously when we assume a sharp drop in the pressure. The second situation where 
sharp transitions are of interest in cosmology are phase transitions like the QCD transition 
or a GUT transition. In these transitions the pressure is continuous, but its derivatives may 
be discontinuous. An example is provided by a first order QCD transition where the sound 
velocity may jump 



Recently, Deruelle and Mukhanov |TT[ derived the joining conditions for scalar met- 
ric perturbations in a spatially flat FLRW model. For perfect fluids (no scalar fields) 
gauge-invariant joining conditions for cosmological perturbations have been derived before 
by Hwang and Vishniac J26j. The difficulty to state the correct joining conditions arises 



because the physical hypersurface of the transition is not necessarily that of constant coor- 
dinate time. In the equality and reheating transitions the physical hypersurface is the one of 
constant density contrast. The gauge-invariant variables describe zero shear hypersurfaces 
as constant time hypersurfaces. Below we derive the joining conditions for general, spatially 
non-flat metric perturbations. Although the method of Deruelle and Mukhanov might be 
simpler than ours, we think that it is worth to view the problem from a different perspective 
below. 

Let us start with exposing the method in general. Assume the spatial transition hy- 
persurface S is defined by its normal n M . In order to join two space-time manifolds along 
S without a surface layer two conditions have to be met ||27||: The induced spatial metric 



hij = gij + riiUj and the extrinsic curvature should be continuous on E. The extrinsic 
curvature is defined as: 

Kij — —£ n hij, (5-1) 

where C n denotes the Lie derivative with respect to the normal n M . In order to compute 
the system of coordinates (i.e. the gauge) and the vector (i.e. the surface of transition) 
have to be specified. Different choices for n M lead to inequivalent junction conditions. Our 



derivation of the joining conditions differs from the derivation in |pT| , where the joining 
conditions are calculated in a coordinate system adapted to the surface of the transition. 
For a more general coordinate system the joining conditions have been obtained by a gauge 
transformation. 

As a simple example we can apply the previous rules to the background model. The 
surface E is defined by qo(rj ) = and the components are given by: n = —a, rii = 0. 
It is then straightforward to show that K z j(rj) = —H{r])5 l j. The continuity of the induced 
spatial metric leads to \im t ^ [a(ri a + e) — a{j] a — e)] = [a}± = whereas the continuity of the 
extrinsic curvature amounts to [a']± = 0. From the Friedmann equations we see that the 
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energy density can not have a jump, whereas the pressure may jump. Let us turn now to 
the case of scalar perturbations. 



A. Scalar perturbations 

The perturbed transition hypersurface is defined by go(^s) + <5<z(?7e,2^) = 0. From the 
last expression we immediately get that the transition now occurs at time rfe = 770 + St] = 
Vo ~ Sq/q'o- I n addition, the normal of £ now reads: 

n o = -a(l + 0), rii = -a^-^- , (5.2) 

% 

where the scalar perturbations are parameterized by the line element Q2.1|) . First we must 
write that the perturbed induced metric hij is continuous on E. Expressing this condition 
for diagonal and off-diagonal terms leads to: 

MVx)]± = iip + H5q/q> } ± (v) = 0, (5-3) 
[E(rm)]± = [E(v)]± = 0. (5.4) 

Second, we must compute the perturbed extrinsic curvature for the vector whose components 
are given in Eq. (|5.2|). We obtain the following result: 



SIC fa) = + H4>)#i + (B-E' + 8q/q' Q )\] , (5.5) 

and therefore, on the surface of transition E, 8K l j takes the value: 

SK'fax) = -{W + (W - H^Wj + {B-E' + 8q/q' )%} . (5.6) 

a q 

Let us notice that the extrinsic curvature on the hypersurface E is a gauge-invariant quantity. 
The two other (gauge-invariant) junction conditions are easily deduced from the previous 
expression and read: 

W + H<P+ (W - H 2 ) 6 4}±(ri) = 0, (5.7) 

Qo 

[B - E> + 5 4) ± ( V ) = 0. (5.8) 
% 



These are the same conditions as obtained in Ref. ||11|| . We have shown that these are also 
valid for all (K = 0, ±1) FLRW models. 

We can also establish what are the junction conditions if one chooses to match the 
perturbations on a surface of constant time as it was done in Ref. | pL Q|| - Since Sq no longer 
depends on spatial coordinates, di5q = 0, the normal to the surface of transition is now 
given by: 

n = — a(l + 0), rii = 0, (5.9) 
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and the extrinsic curvature takes on the form: 



SK^irj) = -[W + Hfi&j + (B- E'f^] . (5.10) 



Requiring the continuity of hij and K l j given by (|5.10D at r/s = rjo, a surface of constant 



time, and using Eq. (|3.3|) leads to the joining conditions in the synchronous gauge: 

[h]± = [h'} ± = N± = [h\] ± = . (5.11) 



These are exactly the conditions that have been used in Ref. [ID|. The joining conditions 



K3|), (|5.4j) , ( |5.7|) , and (|5.8| ) are not equivalent to the conditions (|5.11| ). This confirms that 
for a sharp transition the choice of the surface of matching is crucial. 

However, all joining hypersurfaces are equivalent if \p]± = 0. This can be easily seen by 
considering a hypersurface given by (|5.2[) in the longitudinal gauge and then make gauge 
transformations to any other gauge. Thus, _E?( LG ) = _E( LG ) = and [5q^ LG ^ / q' ]± = from 
( j5.8|) . From [p]± = and the Friedman equations [H']± = follows. The joining conditions 
(PP and (|3) imply 

[^,( L G)] ± = [^,(LG)/ + W 0(LG)] ± = (5 _ 12) 

For a perfect fluid the equation of motion </>( LG ) = ^( LG ) reduces the joining conditions to the 
continuity of the metric perturbations and its derivatives. Consider now all hypersurfaces 
that are related to the zero shear (longitudinal gauge) hypersurface by the gauge transforma- 
tions (2.2) where G C^. Then from the gauge transformations ( |2.2|) and the junction 



conditions it follows that all metric perturbations and its derivatives have to be continuous 
in any gauge that is smoothly connected to the longitudinal gauge. Therefore, all these 
hypersurfaces are equivalent to the constant time hypersurface (|5.9| ) if the pressure does not 
jump. We conclude that Grishchuk's joining conditions are correct, provided [p]± = 0. We 
argue in the next section that this is the case in his reheating transition. 



B. Vector perturbations 



Let us define the line element for vector perturbations to read (here again we use the 
notations of Ref. [0 except for the signature of the metric): 



ds 2 = a(r]) 2 {-drf - 2S , i dz i d77 + by + + F^dx^} , 



(5.13) 



with Si and Fj being transverse vectors, i.e. S { 



0. Fictitious perturbations can be 



generated by performing the infinitesimal change of coordinates: 

77 = 77, x i = x i + C i , (5-14) 
where £\ = 0. Under this transformation Si and Fi change according to the equations: 

Si = Si-Ci, F i = F i + ( i . (5.15) 
Therefore, we introduce the gauge-invariant dragging potential ]TJ| defined by 
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= Si + F/ 



(5.16) 



Since there is no possible vector contribution to the normal of the spatial hypersurface 
E, the components of the normal are simply those of the background, i.e. no = —a, n, = 0. 
From the continuity of the induced metric we get that [Fi]±(r]) = and from the expression 
of the extrinsic curvature 

8K i J = ±(E j * + E< {j ), (5.17) 

we deduce that the second junction condition is: [Hj] ± (r;) = 0. These junction conditions are 
gauge-invariant. In the synchronous gauge, they simply reduce to the continuity of the metric 
and its derivative. In the gauge-invariant formulation we have to demand the continuity of 
Hj only. The fact that there is only a single condition reflects the very different behavior 
of rotational perturbations compared to density perturbations or gravitational waves. For 
perfect fluids the equation of motion for H, reads 

Z' l + 2HZ l = 0, (5.18) 

giving rise to a decaying solution oc I /a 2 fixed by one initial condition only. 



C. Tensor perturbations 

Finally, we treat the case of the gravitational waves. The perturbed tensor line element 
is given by ( j2.15|) . As for the vector case, the normal to the surface of transition is simply 
the one of the background. The continuity of the induced metric leads to [hij]±(rf) = and 
a straightforward calculation of the extrinsic curvature, 

5K^ = (5.19) 

shows that the derivative of the metric must be continuous as well, namely [/iy]±(^) = 0. 
There are no gauge dependences in this sector anyhow. 



VI. A SMOOTH TRANSITION 



In this section, we analyze the approach taken by Grishchuk in Ref. [fT0[] . He assumed 
that during inflation, the scale factor is given by a oc |?7| 1+/3 with 1 + (3 < 0. During this 
stage, the function 7(77) is constant and equal to (2 + /3)/(l + (3). Then, instead of matching 
directly inflation to the radiation stage characterized by a(rj) oc (r]—r) e ), 7 = 2, he introduced 
a smooth transition in between. Physically, this smooth transition represents the reheating 
of the Universe. It begins at 77 = rji — e (the end of inflation) and ends at 77 = r/i + e (the 
beginning of radiation). Note that the parameter e introduced above is different from the 
one used in Ref. Here, e is small compared to 771 because we assume that reheating is 
fast. In the limit e goes to zero, we recover the sharp transition considered before for which 
7(77) becomes an Heaviside function jumping from (2 + /3)/(l + (3) to 2. In the case of a 
smooth transition, without taking into account all the details of the reheating process, we 
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do not know how the scale factor evolves between 771 — e and 771 + e. It is clear that the 
function 0(77) [and therefore 7(77)] is probably complicated for this stage of the evolution. 
The idea of Ref. |L0] was to assume that the function 7(77) is given by: 

7 („) = i±^_ + _£_ tanh , (6.!) 



2(1+13) 2(1+13) 

where s is a parameter controlling the sharpness of the transition. This equation holds for 
inflation and reheating, i.e. for rj between —00 and 771 + e. Such a postulated behavior 
for 7(77) can be justified with the help of the following two arguments. First, we have 
7 ~ (2 + /3)/(l + 0) when 77 < 771 — e (the parameter s must be chosen such that the 
tanh reaches quickly the value —1; it is sufficient to take |e/s| ^> 1) and we recover the 
fact that 7(77) is constant during inflation. During reheating 7(77) smoothly passes from 
(2 + /3)/(l + (3) to its exiting value 7(771 + e) ~ 2. Second, more physically, 7(77) is related to 
w by the equation p/p = w = — 1 + (2/3)7(77). The introduction of the expression ( |6.1| ) in the 
last equation reproduces the expected behavior of to in a reasonable model. Therefore, Eq. 
flS.lp gives a reasonable approximation of the real (exact) complicated function 7(77) even if 
details of the reheating process cannot be taken into account in such a simple approach. 

From the previous discussion, it is clear that 7(77) is always a continuous function, even 
at 77 = 771 + e where the explicit joining was performed in Ref. IlIJQ- This means that [p]± 



vanishes, see Eq. ( |3.1|) . In Ref. ||11|| , Deruelle and Mukhanov criticized the calculations done 



in Ref. [10] by means of the smooth transition described before, arguing that the junction 
conditions were not taken into account properly. We have shown in the previous section 
that if the pressure is continuous, the two sets of matching conditions, Eqs. (|5.3|), fl5.4|), 



1|) and Eqs. ( |5.11 ) are equivalent. Therefore, the claim of Deruelle and Mukhanov 



is not appropriate. For a smooth transition, the matching conditions used by Grishchuk are 



perfectly justified since they coincide with the ones derived in Ref. Jll|. The argument of 
Deruelle and Mukhanov would be correct if the transition were sharp and 7(77) discontinuous 
at 77 = 77! + e. Moreover, the exiting values of the functions 7(77) and 7' (77) at the joining 
reheating-radiation are: 

7(77! + e) = 2, 7 , (77 1 + e) = 0, (6.2) 



in contradiction with the claims of Ref. but in accordance with what is written by 
Grishchuk. 

The next step would be to solve Eq. (|3.24|) for 7(77) given by Eq. (|6.1| ). It was shown in 
Ref. ]TD| that the integration of Eq. ( |6.l| ) can be performed and provides us with the function 
a(j]). However, obtaining the corresponding 0(77) is not possible. This is not a problem since 
the potential (a^/7)" '/ \a^y) depends only on 7, a and their derivatives. However, even the 
simple form (|6.1|) is to complicated to allow a direct integration of Eq. ( |3.24| ). Nevertheless, 



we can follow the evolution of p through inflation and reheating. For 77 < 77! — e, 7(77) is a 
constant and the equation ( |3.24| ) can be solved. The solution reads: 



In this paper 771 was used instead of 771 + e to denote the end of reheating. It is very important 
to distinguish these two events. 
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ii = {kn^lAxJp+^kn) + A 2 j_^{kn)] ■ (6.3) 

This solution is the same as for gravitational waves. The initial conditions A 12 are fixed 
by the quantum-mechanical generation of the density and metric fluctuations: see the Ap- 
pendix. The Ai t2 differ from the A\$ and A\p introduced previously. From Eq. (|6.3j ), we 
can determine the value of /x(t/) just before reheating: 

aim -e)~ - l Al [k(rn - e)]^ 1 ~ a 1 Al (km) p+1 , (6.4) 

m/ ' 2^r(/3 + |) LW n 2^r(/5 + |) 1 /; 1 ; 

because krji <C 1 and e <C rji. Between rji — e and T]i + e the function 7(77) is no longer 
a constant and the solution ( |6.3|) can no longer be used. In order to evolve /1 through the 
reheating transition we use the superhorizon solution \i ~ a^/j to obtain 



M(m + e) - — = =«fa + ^)\flijh + e) - Kvi - <0+/— > ( 6 - 5 ) 

«(m - e) ^7(^1 - e) V 7* 

because a(r/i + e) ~ a(?yi — e). 7, is the value of 7(77) during inflation. This relation should be 
compared to Eq. (81) and to the relation (J,\ m -o = A 4 1 771+0 below Eq. (48) of Ref. fli~0| . From 



Eq. ( |6.5|) ) it is clear that the ratio ^(771 + e)/fi(r]i — e) is not 1 but proportional to l/y/ji 
This factor is huge when 7^ is close to (de Sitter case). Therefore, the mistake in Ref. flt| 
was not due to the use of wrong junction conditions but to the fact that the function /i(r/) 
was not evolved correctly through the reheating transition: actually 7(7/1 — e) 7^ 7(7/1 + e) 
implies //(771 — e) ^ u(r)x + e). 

The value of $ for superhorizon modes during inflation is obtained by inserting ( |6.3| ) 
into ( |3.27| ), using the Taylor expansion of the Bessel functions and substituting (|6.4j) in the 
corresponding expression. The result reads: 

$, * -4±l^ ^- £) . (6.6) 
2/3 + 3 2 a(r/i) V ; 

This equation will be used below. 

Let us turn now to the second transition, i.e. the transition radiation-matter taking place 
at equality. The matter era is described by 0(77) oc (77 — 7/ m ) 2 and 7 = 3/2. In principle, one 
should do the same interpolation at equality as was done at reheating. In Ref. JTIJ, this was 



not done and the second transition was treated as a sharp transition. As was shown in |jTT| . 
the joining conditions (|5.11|) are not correct in general (i.e., for any residual gauge fixing). 
However, if one specifies the synchronous gauge in the matter dominated epoch to be the 
comoving one, then the joining conditions ( |5.11| ) are fine at the equality transition |TT] . 



That is because the density contrast vanishes at the leading order [i.e., it is proportional to 
(krj) 2 ]. This was actually done in [10|. Neglecting the decaying mode leads to the solution 



h = C 1 , / i / = T ^CiA; 2 (7/-7/ m ) 2 . (6.7) 

The coefficient G\ is related to \i{r\\ + e), the value of \i at the end of reheating, by (see Eq. 
(84) of Ref. |3): 
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Therefore, everything is known at the matter stage. Let us calculate the Bardeen's gauge- 
invariant potential. Inserting Eq. ( |6.7| ) in the formulae ( |3.13| ) - ( |3.15| ) gives: 



'To Cl 



-//(?7i + e). 



10V2a(ryi)' 

We may now take into account Eqs. ( |6.6|) and (|6.9| ) to arrive at 



(6.9) 



$7 



2(3 + 3 



2 2/3 + 3 1 

_ rs^ 

1 + /3 fi{vi -e)~5(3+ll + Wi 



(6.10) 



Thus, we obtained the 'standard result' ( |4.7| ) entirely within the synchronous gauge, without 
any reference to the constancy of ( or the joining conditions of Deruelle and Mukhanov. From 
the last expression, it is clear that it is crucial to evaluate correctly the ratio /i(?7i + e)/ ji{rji — 



Finally, it is also interesting to calculate the ratio of $ and /i gw at superhorizon scales 
today. This quantity is of relevance for observations since it is related to the ratio of grav- 
itational waves to density perturbations contributing to the CMBR quadrupole anisotropy. 
The leading term of the 'growing' mode of the gravitational waves is constant in time, thus 
its value today is equal to its value at time r) — rji — e. Therefore we obtain: 



h gw (today) 



A? tin 



a(r/i - e) A 1 a(r/i) 



(6.11) 



where we used the fact that the equations of motion of /i, Eq. ( |3.24| ), and /igw, Eq. ( |2.16| ), 
are the same during inflation (as long as 7 is almost constant) until the onset of reheating 
at 7]i — e. From Eqs. ( |A24| ) and ( |A26| ) of the Appendix, it follows that the rms amplitude 
of gravitational waves, which takes into account both polarizations, reads: 



h ~ 

"rms — 



V2 

71 



tiVi ~ e) 



&2 



On the other hand, Eqs. ( |33 



^rms — 



H) and (TAT5D imply that: 
3 



(6.12) 



[J 


tivi - e ) 


1 li 





(6.13) 



Therefore, the value of the ratio /i rms /$ r ms today is: 
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At 



\/6 



Wi 



(6.14) 



where the equations 



of the Appendix have been used. This finally proofs that the 



closer the inflationary epoch is to the de Sitter space-time, the less important are large-scale 
gravitational waves in the CMBR today. 
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VII. CONCLUSION 



In Ref. [[TIJ Grishchuk claimed that the magnitude of superhorizon scalar metric pertur- 
bations is most likely smaller than the amplitude of superhorizon gravitational waves. We 
have shown in Sec. VI that this result is wrong, because the time evolution of the scalar 
metric perturbation through the (smooth) reheating transition was not calculated correctly. 
With the appropriate correction we recover the 'standard result' for the rms amplitudes at 
superhorizon scales 



h r 



20 f— m P iV v 



(7-1) 

slow— roll 



today V / 6 V 1 V 

where the slow-roll approximation is valid if 7$ <C 1. However, in the limit 7, — > linear 
perturbation theory breaks down since Eq. ( |6.13| ) blows up. Thus, for power-law inflation, 
the model which we considered in detail, one cannot make the slow-roll approximation 
arbitrarily precise by making 73 arbitrarily small. 

Recently, Deruelle and Mukhanov corrected the result of |10| within the framework 
of sharp transitions. We rederived their joining conditions in Sec. V and extended them 
to non-flat FLRW models (K 7^ 0). Moreover, we derive the joining conditions for the 
vector and tensor perturbations. According to Deruelle and Mukhanov, Grishchuk made 
two mistakes: He took the wrong joining conditions and he used the wrong equation of state 
(expressed in terms of 7) at the reheating transition. However, Grishchuk introduced a tanh 
to interpolate the pressure between the inflation and radiation epochs. Therefore, both his 
joining conditions at the reheating transition and the equation of state after reheating have 
been used correctly. 

A commonly used derivation of the 'standard result' has been criticized by Grishchuk in 
Refs. [|T^,|rjJ]. This derivation is based on the conservation of certain quantities for superhori- 
zon modes. These 'conservation laws' are essential in Refs. []6|,p2|-f24| . We have investigated 



Grishchuk's arguments in Sec. Ill and IV, where we have shown that his criticisms are not 
correct. 

Note added: After our work was finished, a paper by M. Goetz (astro-ph/970427) ap- 
peared. In this paper he independently reaches one of the conclusions of our paper, namely 
that Grishchuk's claim on the emptiness of the 'conservation law' is wrong. 
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APPENDIX: QUANTIZATION 

In this appendix, we briefly review how the perturbations are generated quantum- 
mechanically in the early Universe. This mechanism fixes the initial conditions, i.e. the 
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coefficients Af, A| w and A\, A 2 in Eqs. ( |2.17 ) and (|6.3|) . It has been emphasized in the 
text how crucial the precise values of these coefficients are to obtain of the final (standard) 
result. 

Let us first consider density perturbations. The normalization of the (perturbed) scalar 
field operator is fixed by the uncertainty principle of Quantum Mechanics. In a spatially 
flat FLRW model, this leads to the following expression: 

S0M - ^ / dkfcfok)^ = ^J^m J ^(v)^ + cUrj)e-^% (Al) 



where Ck and c k are the annihilation and creation operators satisfying the usual commutation 
relation. This equation agrees with Eq. (95) of Ref. [10]. The initial conditions are deter- 
mined by demanding that at some initial time rjo (e.g., at the beginning of inflation), the 
scalar field be placed in the vacuum state: Ck(r/ = r)o)\0) = 0. Then, solving the Heisenberg 
equation of motion shows that the operator (pi (77, k) can be written as: 

- / i\ r , ,U k + vt + .U* k + V kl ... 

^"■ k) = ^ W "» ) W + cLkM W 1 ' (A2) 

where the function (uu + v^)(t]) satisfies: 

(u* + v* k )" + (k 2 - -)(u k +v* k ) = 0. (A3) 



The initial conditions translate into the statement that Uk{i]o) = 1 and Vk{i]o) = 0. In the 
high frequency regime, this implies the asymptotic behavior lirm : _ ) .-|- 0O ("Ufc +v£) = e~ tk(j) ~ m \ 
The normalization of the perturbed scalar field fixes automatically the normalization of 
the scalar perturbations of the metric since they are linked through Einstein's equations. In 
the high frequency limit, this link is expressed through the formula: 

lim jl(r), k) = — v / 2Ka<^i(?7, k) , (A4) 

which follows most easily from Eqs. (|3.26 ) and (|3.25|) . This equation is the same as written 

at the bottom of p. 7168 of Ref. ||10|| . This allows to find immediately the asymptotic 
behavior of the operator /l(ry, k): 

e — ifc(?7— ?yo) gife(7?-J7o) 

lim #(77, k) = -4v/^pi[c k (r/o) j^— + cL k M — /=H, (A5) 



where lp\ = (Gh) 1 ^ 2 is the Planck length. For simple models where the scale factor is 
a (v) = ^oM 1+/3 ; the exact solution for fi(rj, k) can be expressed as: 



/}(r/,k) = (A;r / )^[A 1 (k)J /3+1/2 (^) + A 2 (k)J_ (/3+1/2) (A ; r / )], (A6) 

where J±(p+i/2) is a Bessel function of order ±(/3 + 1/2). The two last equations imply that 
the (so far arbitrary) operators -4i(k) and A 2 (k) are given by: 
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iy/8irlpi 1 



ii(k) = i^^ iikV ° +7f)c ^o) ~ e~^ + W k M], (A7) 



i 2 (k) = -^^^=[e^-^c k (Vo) + e-«**-*V_M]. (A8) 

COS p7T a/2k 



Essentially, these relations agree with Eqs. (102) of Ref. JlC . 

The power spectrum for the Bardeen potential can now be computed. Its definition is 
given in terms of the two-point correlation function for $(77, x): 

(0|$(r /) x)|.(r /) x + r)|0)= f° ^ S -^-k s P^(k). (A9) 

jo k kr 

Using Eq. (3.27) which expresses the link between \x and $, we find the following expression 
valid for the simple models evoked previously and for long wavelengths: 

P * {k) = f2wJp^ + 5/2) (A10) 

So far, this is the result during inflation on superhorizon scales. Although no explicit time 
dependence is visible in ( |A10| ), the value of Pq>{k) changes during reheating and during 



the equality transition. In particular, the spectrum today is equal to: P$ (today, k) 
T(k)P<s> (initial, k). For super horizon modes, the transfer function is given by: lim^o T(k) = 
[3(2/5 + 3)] 2 /[5(l + P)^] 2 . If space-time during inflation was close to a de Sitter phase, then 
the spectrum today is the well-known Harrison- Zeldovich (scale-invariant) spectrum (i.e. 
13 < -2). 

This result permits us to carry the quantum-mechanical initial conditions to the classical 
level. Let us define the 'classical spectrum' for the classical quantity $(77, x) by the following 
expression: 

1^1 \*.t w Iv dx$(r/, x)$(r/, x + r) /»,,<, 
$ 77, x)$(t/, x + r = ^ W, r ' Kh '- All 



where P${k) is given by: 



^(*) = Al*fa,k)| 2 - (A13) 

Z7T Z 



In this expression, $(77, k) can be calculated using Eqs. (3.27) and (6.3) where the un- 
known coefficients A 1: A 2 appear. Requiring that P$(k) = P$(k) fixes This last 
equation relies on the ergodic assumption, namely that ensemble averages are equal to 
spatial averages. In addition, if we have the following behavior for the classical ^(77, k): 
lim fc ^ +00 /i = — 4^lpie~ tk ( v ~ m ^ / V2k as it is suggested by Eq. (|A~5|), then the classical 
initial conditions A±, A2 are completely determined. They read: 

Al = ZV8 7 16 , A 2 = iA x e~^. (A14) 

cos p7i V2/c 



25 



Finally, we define the root mean square value $ rms by: 



$rms = ^k 3 P„(k). (A15) 

This quantity is used at the end of Section VI. 

Let us now consider gravitational waves. The problem is very similar to the previous one 
since the quantization of gravitational waves is equivalent to the quantization of two scalar 
fields (representing the two independent degrees of freedom of the wave). Assuming again 
that the initial state is the vacuum leads to the following equation for the gravitational wave 
operator: 

kiiV, x) = j^j~ 2 £ / dkp^mUv, k K k ' X (M6) 

+ « + 4K t (^)e- ik - x ]. (A17) 

In these formulas, pfj(k) is the (transverse-traceless) polarization tensor and the summation 
over s represents the summation over the two states of polarization of the wave. The 
polarization tensor is normalized as pfj(k)p s u (k) = 25 ss . The function u s k + vf.* satisfies Eq. 
(2.16) (and in fact does not depend on the state of polarization s). The initial conditions 
are: u s k (rjo) = 1 and v s k {rjo) = 0. This implies the following asymptotic behavior for the 
operator fi s gw (r],k): 



p -ik(r)-r) ) p ik(r)-r) ) 

lim p,' fa k) = V^PiKfo) + 0^)—=-]. (A18) 

fc ^°° V2/c \/2k 



This equation is similar to Eq. ( |A5|) . Since the exact solution for ft s gw (r], k) can be expressed 
as: 

^ w (r?,k) = (k V )V 2 [ArOL,s)J p+1/2 (k V ) +Ar(k,s)J^ +1/2) (k V )} , (A19) 
this fixes the operators Af w (k, s) and A| w (k, s): 

if(k, s) = -^^^=[e^ + ^c^) - e'^^M (A20) 

COS LsIT \f 2i\\i 



if(k, s) = ^ * [e^-^fa.) + • (A21) 



The result for gravitational waves is very similar to the result for density perturbations. 

We can also compute the power spectrum for gravitational waves (it involves the calcu- 
lation of (0 1 /i^j- (?7, x)/i 4J (r7, x + r)|0)). For long wavelengths, we obtain: 

I 2 1 

Ph(k) = f 2^-2 cos 2 ( /?vr )P(/5 + 3/2) (A22) 

Interestingly enough, the /c-dependence of the density perturbations and gravitational waves 
power spectrum is the same. This behavior is a special feature of power-law inflation. As 
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for the density perturbations case, the calculation of the classical spectrum leads to the 
determination of the classical initial conditions. The classical spectrum is given by: 



(/lyfox^'fox + r)) 
where P^{k) can be expressed as: 



00 dk sin kr 
k kr 



/4w0/; k ) 



(A23) 



(A24) 



Requiring Ph{k) = Ph(k) and demanding the same high frequency behavior for the quantum 
and the classical solution leads to the determination of Af w , Af which appear as free 
numbers in the expression of /i gw (?7, k), Eq. (|2.17|) . We obtain the following result: 



gwi 



l^ll 



\AT\ = \A 2 



(A25) 



Thus we have derived the result used in Sec. VI: |Af w /Ai| = 1. 
Finally, the classical quantity h rms is defined by: 



h r 



k 3 P h (k). 



(A26) 
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TABLES 



H = a' I a 
&B,il>,E 
$, ¥ 
a 

Si, Fi 
ha 



H = aH, where H is the expansion rate. 
Perturbed scalar metric coefficients, see (|2.1| ). 
Gauge-invariant scalar metric-perturbation, see ( |2.3| ). 
Proportional to <£, see ( 2.12| ). Denoted u in fp~5p| . 
Gauge-invariant velocity potential, see ( |3.25| ), Denoted v in [15] 
Perturbed vector metric coefficients, see ( ^.13[ ). 
Gauge-invariant vector metric-perturbation, see (5.16). 
Perturbed tensor metric coefficient, see (2.15). 



TABLE I. Notation of gauge-invariant formalism. 



a = a 1 1 a a = aH, where H is the expansion rate. 

7 = 1 — a' /a 2 7 is related to w = Po/po by 7 = [3(1 + w)/2](l — K/a 2 ). 

h, hi Perturbed scalar metric coefficients, see (|3.2|) . 

u,v Residual-gauge-invariant scalar metric-perturbation, see ( p. 13 ). 

H Proportional to u, see (HI)- 

TABLE II. Notation of synchronous gauge. 
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